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We examine the notion of nonclassicality in terms of quasiprobability distributions. In particular, 
we do not only ask if a specific quasiprobability can be interpreted as a classical probability density, 
but require that characteristic features of classical electrodynamics are resembled. We show that 
the only quasiprobabilities which correctly describe the superposition principle of classical electro- 
magnetic fields are the s-parameterized quasiprobabilities. Furthermore, the Glauber-Sudarshan P 
function is the only quantum-mechanical quasiprobability which is transformed at a classical atten- 
uator in the same way as a classical probability distribution. This result strengthens the definition 
of nonclassicality in terms of the P function, in contrast to possible definitions in terms of other 
quasiprobabilities. 
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I. INTRODUCTION 

Since the early days of quantum optics, many effects 
have been observed which are incompatible with the clas- 
sical electromagnetic theory of light. Already Einstein's 
explanation of the photoelectric effect [I] , for which the 
photon was postulated, was in contradiction to the clas- 
sical theory and gave rise to quantum mechanics. Much 
work has been done to exactly describe the difference be- 
tween the classical electromagnetic field and the quantum 
optical counterpart. 

In the 1960s, physicists began to treat this problem in 
terms of phase space distributions. Classically, any state 
of a single mode of light can be described by a classical 
probability distribution over complex field amplitudes a, 
which span the so-called phase space. Quantum mechan- 
ically, phase-space distributions are not uniquely defined. 
Well-known examples are the Wigner distribution [5] , the 
Glauber-Sudarshan P function [3j|4], and the Husimi Q 
function [5]. They have been generalized as members 
of the class of s-parameterized quasiprobability distribu- 
tions [5] . The most general form of quasiprobabilities has 
been given by Agarwal and Wolf [7]. 

Nonclassicality can be defined by comparing such a 
quasiprobability with a classical probability distribution. 
However, this definition depends on the particular choice 
of quasiprobability. The generally accepted definition of 
nonclassicality is based on the Glauber-Sudarshan dis- 
tribution P function [HI [S]. This definition is designed 
in such a way that any correlation function of a classical 
state can be modelled by classical electrodynamics. How- 
ever, the P function may be highly singular, such that it 
cannot be directly reconstructed from experimental data. 
Therefore, many nonclassicality criteria had to be devel- 
oped. The most general, but still simple criterion is based 
on so-called nonclassicality quasiprobabilities [10] , which 
has already been applied in practice [TTJ [T2]. Another 
frequently used sufficient criterion is based on the nega- 
tivity of the Wigner function |13j . 

In this article, we strengthen the role of the P func- 
tion as the key quantity for the definition of nonclas- 



sicality. The foundation of our considerations is the su- 
perposition principle of the classical electromagnetic field 
of light, which is employed at every beam splitter in an 
optical experiment. We demonstrate that not every gen- 
eral quasiprobability is transformed by a beam splitter 
in the same way as a classical phase-space distribution. 
More precisely, only the s-parameterized quasiprobabili- 
ties behave classically in this situation. Hence they are 
the only candidates for a reasonable definition of nonclas- 
sicality. Furthermore, we show that only the Glauber- 
Sudarshan P function behaves under the action of an at- 
tenuator like a classical phase-space distribution. Hence, 
it is the only quasiprobability which completely resembles 
the classical behavior of linear operations on the classi- 
cal field of light. This shows that defining nonclassical- 
ity in terms other quasiprobability distributions, e.g, the 
Wigner function [14], leads to contradictions in the in- 
terpretation of the classical character of an experiment. 

The paper is structured as follows: In Sec. II, we com- 
pare the action of a beam splitter on classical phase 
space distributions and quantum-mechanical quasiprob- 
abilities. In Sec. Ill, the attenuator is analyzed for clas- 
sical and quantum optics. In Sec. IV, we give an explicit 
example to demonstrate that only the P function can dis- 
tinguish between classical and nonclassical behavior cor- 
rectly, and show that a state with a nonnegative Wigner 
function may lead to measurement results which cannot 
be explained by classical electrodynamics. 



II. ACTION OF A BEAM SPLITTER 
A. The classical superposition principle 

Let us consider the electric field of a monochromatic 
light wave with a fixed frequency uj. The spatial depen- 
dence of the field shall be neglected in this manuscript. 
Then, the field can be fully characterized by a complex 
amplitude a 

E{t-o) = ae lulT + a* e-^ T . (1) 
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FIG. 1. Scheme of a beam splitter. 



The maximum amplitude is given by E m&x (a) — 2Re(o), 
and the phase can be specified by arg(o). The set of all 
complex o is referred to as phase space. 

A characteristic feature of light is the possibility of 
linear superposition of electromagnetic fields, which can 
be achieved by a beam splitter, see Fig. [T] The relation 
of input and output fields E(t, ai) are given by 



E(t, a 3 ) 
E(t, o 4 ) 



U 



E(r,ai) 
E(t, a 2 ) 



where the matrix 



U = 



t* 



(2) 



(3) 



has to be a unitary one. Throughout the manuscript, we 
set the global phase ipu = 0. Instead of considering the 
fields, this map can also be considered in terms of the 
complex amplitudes Oj, 



(14 



U 



a 2 



(4) 



If all parameters are non-zero, the knowledge of two arbi- 
trary amplitudes o^ , aj complete determines the missing 
two amplitudes. This is trivial if the amplitude oi, a 2 or 
03,014 are given. If 01,03 are known, the missing two 
fields are determined by. 



Q'2 
O4 



-t 



Oil 
0:i 



All other cases can be treated analogously. Therefore, 
whenever two classical fields E(t, o^), E(t, aj) are given, 
there exist two additional classical fields such that Eq. ([2| 
is satisfied. 

Moreover, we may consider a statistical ensemble of 
fields, modelled by a probability distribution P c i.(o). In 
this situation, every physical quantity F(a) has to be 
averaged over this probability distribution, 



(F) 



d 2 oP cl .(o)F(o) 



(5) 



Furthermore, the action of a beam splitter on an two- 
mode input field characterized by its joint probability 
density P12 (01,02) is given by 



Therefore, for each classical input field ensemble there 
exists a classical output field, defined by the probabil- 
ity distribution P34 (03,04). Due to the linearity of the 
relations, this holds for an arbitrary choice of two given 
modes i,j. 



B. The beam splitter and quantum phase-space 
functions 

Quantum mechanically, the expectation value of an ar- 
bitrary operator F can be calculated in the form [TS] 



(F) = / d 2 aP n {a)F h {a). 



(7) 



Here, the so-called quasiprobability Pn{a) is given 
by the Fourier transform of its characteristic function 
<M/3) (cf. 0), 



Pn(a) = 



1 



d 2 /3$ n (/3)e 



(8) 



which itself is directly connected to the quantum state of 
the oscillator, described by its density operator p: 



$o(/3)=Tr(pe' 3at -' 3 * a )fi(/3,/3*). 



(9) 



The filter function Vt(f5) characterizes the special kind of 
the quasi-probability, which we examine more in detail 
below. The function F^(a) is connected to the operator 
F in a similar way, with f2(o) being a suitable inverse 
filter. 

The different known types of quasiprobabilities are de- 
termined by a specific filter f2(/3). For instance, the 
Wigner function can be obtained by fJ(/3) = 1, the P 
function by fl(l3) = e'@' / 2 and the Husimi Q function by 
f2(/3) = e~'^' I" 1 . The latter quasiprobabilities are gener- 
alized as s-parameterized quasiprobabilities in [B], where 
the filter has the form f2(/3) = e s '^l / 2 . In general, we re- 
quire that S7(/3) is an analytic function in (3 and j3* , has 
no zeros, and satisfies O(0) = 1. Under this conditions, 
we derive our first important result: 

Theorem: Only the s-parameterized quasiprobabili- 
ites resemble the classical action of a beam splitter. 

Proof: It appears that the action of the beamsplitter 
on ensembles of electromagnetic fields can be examined in 
a convenient way in terms of the characteristic function 
'^(/^j ^4) of the probability distribution which is 
defined as 

*34(/33,/3 4 ) = J d 2 a 3 J d 2 o 4 P 3 4(a3,a4) 

Then Eq. @ can be rewritten in terms of characteristic 
functions as 



P3 4 (o 3 ,o 4 ) = Pi 2 (t*a 3 - ro4,r*o 3 + to 4 ). (6) 



$34(/33,/3 4 ) = <Pi 2 (t*p 3 -r/3 4 ,r*p 3 + tp 4 ). (11) 
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In the following, we seek for the characteristic functions 
^n(^3: Pa) °f the quasiprobabilities which also satisfy 
this equation. The theorem states that only quasiproba- 
bilities with n(/3) = e s '^' I" 1 fulfill this requirement. 

To show this, we first note that in the case of two-mode 
fields, the filter function of the two-mode characteristic 
function ^n(^3i fit) has to appear as a product of single 
mode filter functions, 

Mft.Ai) = Ti(pe^+^-^-^)n{0 3 )n(^). 

(12) 

Then, the action of the beamsplitter on the creation 
and annihilation operators can be written analogously 
to Eq. @, cf. [H]: 



fl.-i 

&4 



u 



«2 



(13) 



Let us denote the characteristic function of the Wigner 
function, defined by f2(/3) = 1, as (As, /?4). We can 
easily show that Eq. (11) holds for $jy (/3 3 , /3 4 ). There- 
fore, it is sufficient to examine the relation of the beam- 
splitter on the filter function fl(/3). More precisely, any 
characteristic function $o(/33, /?4) satisfies the classical 
relation Eq. (11) if and only if the filter function £l((3) 
satisfies 

OCAOfW - - r/3 4 )n(r*/3 3 + tfa). (14) 

The only conditions which are imposed on the beamsplit- 
ter coefficients t, r arise from the unitarity of the matrix 
U in Eq. ^: 

|t| 2 + |r| 2 = l. (15) 

All filter functions f2(a) which are entire in a, a* and 
have no zeros, can be written in the form (cf. [7]) 



rt(P) = exp 



k,l=0 



(16) 



We may insert this expression into Eq. ( 14 ) and compare 



the coefficients of the series in the exponential functions. 

Since we require fi(0) = 1, we need cqq — 0. Let us now 
examine all powers of (3 k {P 3 ) with k + l > 0. Comparing 
both sides of the equation, we find the condition 



c H = c kl W) k t l + {r*ft 



(17) 



The equation is fulfilled for a fixed pair fc, I if the expres- 
sion in brackets equals to 1 for all allowed parameters 
r, i, or if cm = 0. It is easy to show that two special cases 
already lead to a large number of vanishing coefficients: 

• If we insert the valid parameters t = r = l/\/2, the 
term in brackets is given by 



[(£*)¥ + (r*)V] = 2/V¥+ l 



(18) 



This expression can only be equal to one if k+l = 2. 
Therefore, it is necessary that Cf-i = for k + I > 2 
as well as cio = cqi = 0. 



• If we insert the valid parameters t = 1/V2 and 
r = and consider the case k + l — 2 only, the 

term in brackets is given by 

[{t*) k t l + (r*) k r 1 } - (1 - (-l) k )/2, (19) 

which equals to 1 if and only if k = I = 1. There- 
fore, we find C20 = C02 = 0. 

Therefore, after setting en = s/2 we find that only filters 
of the form 



nC8) = exp{ s |/3| 2 /2} 



(20) 

are candidates to satisfy Eq. |l~4"| ). 

As a last step, we still have to show that the filter 
function pOl fulfills Eq. (fTil). To this end, we start with 



the right hand side of the latter equation: 

= exp{ S [(|i| 2 + |r| 2 )|/? 3 | 2 + (|r| 2 + |i| 2 )|/3 4 | 2 

+ (rt - trffapt + (rt ~ tr)f3* 3 ^}/2} . (21) 

Now, we can insert the relations |T5"1 ) and easily verify 
that 



exp { i 



l/3 4 | 2 ]/2} 



(22) 



The last term exactly equals to the left hand side of the 
desired equality, cf. Eq. (14). This concludes the proof. 
□ 



III. ATTENUATION OF THE 
ELECTROMAGNETIC FIELD 

So far, we have shown that the action of a beam 
splitter on quantum mechanical quasiprobability distri- 
butions can only resemble the classical behavior if the 
quasiprobabilities belong to a special class, namely the 
s-parameterized ones. Hence, as the Q-function, Wigner 
function and P function all belong to this class, they are 
suitable for a classical description of the beam splitter 
action. However, there are s-parameterized quasiprob- 
ability distributions which are still not suitable for the 
definition of nonclassicality. For instance, the Husimi Q 
function always satisfies the requirements for a classical 
probability density. 

In this section, we will impose an additional require- 
ment. Classically, damping of the electromagnetic field 
can be considered as the action on the beam splitter, 
where a signal field ai is overlapped with the vacuum 
field ai — 0, leading to a resulting attenuated field 
«3 = tax. The parameter 77 = \t\ 2 is commonly referred 
to as (quantum) efficiency, which diminishes the ampli- 
tude of the signal field. The joint probability density of 
the input field is given by the product of the probability 
densities of the state of interest and the vacuum field. 



Pi 2 (ai 1 a 2 ) = Pi{ai)5(a 2 ). 



(23) 



The probability density of the attenuated output field 013 
is then obtained as the marginal of ^34(013, 1x4), given by 
Eq. ©: 



^3(^3)— J d 2 a^Pi{t*a^ — ra^)8(r*a^ + ta^ 



= pp P l(«3/*) 



(24) 



Vice versa, given an classical output field distribution 
Pzi&s) an d the quantum efficiency 77 = \t\ 2 > 0, there 
exists a classical distribution Pi(ai) = \t\ 2 P-j,{a\t) for a 
classical input field of the attenuation process. 

Quantum-mechanically, attenuation is also considered 
as overlapping the signal field with vacuum at the beam 
splitter, described by quasiprobabilities Pn.i{ a i) an d 
Pn,vac(oi2)- Therefore, the input state is described by 
the quasiprobability 



Pn, 12(01, "2) = Pn,i{ai)Pn,va.c(a2)- 



(25) 



This leads us to the following result: 

Theorem: The Glauber-Sudarshan P quasiprobabil- 
ity is the only quasiprobability which allows to describe 
the action of an attenuator in the same way as the clas- 



sical electromagnetic theory, see Eq. ( 24 ) . 

Proof: Again, we examine the quasiprobabilities in 
terms of their characteristic functions. The classical re- 



lation ( 24 1 can be rewritten as 



*3(Ai) = 3i (**&). (26) 
On the other side, we have to insert the Fourier trans- 



form of Eq. (25 1 into the quantum-mechanical version of 
Eq. jllj: 

$fi,34(/33,/3 4 ) = *n, 1 (**^3-r J 9 4 )*n,vac(r*j93+«j94). (27) 

Since we consider only the output mode 3, we have to 
take the marginal with respect to the mode 4. In Fourier 
space, this is simply achieved by setting ^4 = 0. There- 
fore, the output mode 3 is described by the characteristic 
function 



*n,3(A0 = *n,i(**/33)*n,™(r*/3 3 )- 



(28) 



Comparing this result with the classical expectation in 
Eq. J26l we find 



4> 



O.vac 



1 



(29) 



as the necessary and sufficient condition for the equality 
of the quantum mechanical and classical quasiprobabili- 
ties. This is only satisfied for the characteristic function 
of the P function, which is obtained by f2(/3) = e'^'"/ 2 . 
□ 

We emphasize that the key argument of the proof is 
that the quasiprobability for the vacuum state exactly 
equals to the classical probability distribution for a field 
with zero amplitude, namely P vac (a) = 8(a). 




FIG. 2. Experimental setup of second-order correlation mea- 
surements. The input field is split by a symmetric beamsplit- 
ter. The intensities of the output field are correlated. 



IV. DISCUSSION 

Let us illustrate our conclusions with an example. 
For this purpose, we consider the second-order corre- 
lation measurement (Fig. ^1), as it is used for photon- 
antibunching experiments |1y| . In classical terms, the 
input field a- ln is split by the symmetric beam splitter 
into two output fields a 1.2 cx whose intensity corre- 



lation is measured, 



G^Vin) 



K| 2 |a 2 | 2 cx|a in | 4 . 



(30) 



Consequently, for any statistical superposition of input 
fields a-,„, the correlation function reads as 



G 



(2) 



Pa.{a)Gf{a)d 2 a 



(31) 



Moreover, one may show that any classical field satisfies 
the inequality [HIIS] 



G { 2 > [G (1) ] 2 , 



with the first order correlation function 



G& 5 = J P c Ua)\a\ 2 d 2 a, 



(32) 



(33) 



describing the intensity of the input field. 

As a first quantum mechanical system, let us consider 
a single photon, which has a Wigner function with neg- 
ativities as well as a highly singular P function. Hence, 
in terms of both quasiprobabilities, the photon can be 
referred to as nonclassical. We obtain the correlation 
functions as 



and 



G (1) cx (fit a) = 1 



G< 2 > cx ((at) 2 a 2 ) = 0. 



(34) 



(35) 



Obviously, the photon docs not satisfy inequality (|32j). 
Therefore, the outcomes of such an experiment cannot 
be explained by means of classical electrodynamics, and 
the term " nonclassicality" is justified. 
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However, the situation changes when we examine a 
slightly different quantum state, namely a single pho- 
ton, which is attenuated by a value of r) = 1/2. Classi- 
cally, the intensity is simply divided by two. Quantum- 
mechanically, the field has to be described by effective 
annihilation operators 

b= (a + a vac )/V2, (36) 



being nonnegative for rj < 1/2. Hence, in terms of the 
Wigner function, one might state that the quantum state 
is classical, although the second-order correlation func- 
tion violates basic laws of classical electrodynamics. As 
we have shown in the previous section, the P function 
is the only quasiprobability which resembles the classical 
behavior of the attenuator, and is free of such a contra- 
diction. 



where a is the annihilation operator of the signal field and 
a V ac is the annihilation operator of an auxiliary vacuum 
field. In this case, the correlation functions are given by 



and 



G (1) oc (S f S) 



G< 2 > cx ((tffb 2 



(37) 



(38) 



Therefore, the classical inequality (32) is still violated 



However, the Wigner function of the input state is given 

by m 



WJa) = -[l-277 + 47 ? H 2 ]e 

7T 



-2|a| 2 



(39) 



V. CONCLUSIONS 

We have examined the action of linear devices, such 
as the beamsplitter and the attenuator, on quantum me- 
chanical phase-space distributions and compared it with 
the situation in classical electrodynamics. Our goal was 
to find out which quasiprobabilities are able to resemble 
the classical behavior, in order to be useful for the defini- 
tion of nonclassicality. We found that for the beamsplit- 
ter all s-parameterized quasiprobabilities, such as the Q 
function, the Wigner function and the P function re- 
flect the classical action. However, the attenuator can 
only be described analogously to the classical picture in 
terms of the P function. Therefore, the latter is the only 
quasiprobability which maps the state of light in the same 
way as classical linear devices do, preserving the classi- 
cal properties of a field. Hence, the P function possesses 
properties which are superior for the definition of non- 
classicality. 
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